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Abstract 

We investigate decay constants of P and D-wave heavy-light mesons 
within the mock-meson approach. Numerical estimates are obtained us- 
ing the relativistic quark model. We also comment on recent calculations 
of heavy-light pseudo-scalar and vector decay constants. 
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I. INTRODUCTION 



Reliable estimates of heavy-light meson decay constants are important, since they 
appear in many processes from which fundamental quantities can be extracted 
Theoretical investigations have focused on estimating decay constants for the weakly 
decaying pseudo-scalar meson and its HQET (heavy quark effective theory) related 
vector meson. Whereas the decay constant of the weakly decaying pseudo-scalar meson 
is of paramount importance for determining fundamental quantities, the decay constant 
of the S'-wave vector meson plays a role in exclusive b — > ulu transitions [0] and in 
radiative leptonic decays of heavy-light mesons . 

While those decay constants have been and continue to be studied intensively, 
the decay constants of the more highly excited heavy-light states have been normally 
ignored. This note attempts to rectify this situation, by predicting decay constants 
for many higher-excited resonances. That could be important phenomenologically on 
several accounts. 

First, CLEO recently observed a significant wrong charm contribution in B decays 

i, 

B(B^DX)^ 10%, (1.1) 

governed essentially by the b — > ccs' quark transitions. [] The B — » DX transitions were 
overlooked in all previous experimental analyses. Under the factorization assumption 
1^, wherein the virtual W ^ cs hadronizes independently to the rest of the system, 
a quantitative modelling of the B DX transitions can be undertaken once theory 
provides the decay constants for D**. 



The prime indicates that the corresponding Cabibbo suppressed mode is included. 
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Second, reliable estimates for decay constants of D** allow one to test whether color- 
allowed and color-suppressed decay amplitudes interfere constructively for the —>■ 
D**^{n^ , , ai , . . .} modes, as has been seen for the D^*^^ {n^ , , ai , . . .} 

transitions 0. Third, such estimates enable us to better predict subtle CP violating 
phenomena. 

Decay constants are defined through matrix elements of vector and pseudo-vector 
currents between meson states and the vacuum. Therefore, in order to calculate them, 
one has to find a way to evaluate hadronic matrix elements. The mock-meson method 
||7HlO| has been frequently used in the literature for that purpose |]7|-[r7|. In this paper 



we follow the same approach, and use the mock-meson method in order to obtain 
expressions for the decay constants of heavy-light mesons, in terms of integrals over 
momentum-space bound-state wave functions. For numerical estimates we decided 



to use the simplest relativistic generalization of the Schrodinger equation |TT|JT8| -pO 
sometimes called the spinless Salpeter equation. 

The rest of the paper is organized as follows. We begin with a brief description 
of the mock-meson method in Section Our approach is based on the j-j coupling 
scheme, since it is more appropriate for heavy-light mesons than the usual L-S scheme. 
Expressions for the decay constants of heavy-light meson states are given in Section 
|T|. The relativistic quark model and our numerical estimates are described in Section 



IV. There we also comment on recent calculations of pseudo-scalar and vector decay 



constants . Our conclusions are summarized in Section 0. 



II. THE MOCK-MESON METHOD 



As already mentioned, the mock-meson approach |]7|-p!0| has been widely used for 
calculations of hadronic matrix elements [[^-p^ . The basic idea of the method is simple. 
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The mock meson is defined as a collection of free quarks weighted with a bound-state 
wave function. The mock-meson matrix elements Ai can then be calculated using 
full Dirac spinors. On the other hand, the physical matrix elements can always 
be expressed in terms of Lorentz covariants with coefficients Ai, which are Lorentz 
scalars. In many simple cases, A4 and A4 will be of the same form. The mock-meson 
prescription then says that in those cases one should simply take Ai = Ai. Indeed, this 
correspondence is exact in the zero-binding limit and in the meson rest frame. Away 
from this limit the mock amplitudes are in general not invariant by terms of order 
pf/ml 

In this paper we are primarily concerned with the decay constants of heavy-light 
qQ mesons. In the m-Q oo limit, heavy quark symmetry tells us that the angular 
momentum of the light degrees of freedom (LDF) in the heavy-light meson decouples 
from the spin of the heavy quark, and both are separately conserved by the strong 



interaction Therefore, total angular momentum j of the LDF is a good quantum 
number. For each j there are two degenerate heavy meson states (J = j ± |), which 
can be labeled as Jj^, where P = (—1)'^^^. This implies that in the case of heavy-light 
mesons the j-j coupling is more appropriate than the L-S coupling scheme. For this 
reason, we first define the LDF states L^) as Clebsch-Gordan (CG) combinations 
of the eigenstates of orbital angular momentum \LMl), and those of the spin of the 
light quark with CG coefficients denoted as C*^^^^ .i^- Combining the LDF states 
with those of the heavy antiquark (with CG coefficients C^^^^j__), we get the qQ 



mock meson state in its rest frame. 
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I -0j-3^^^^r.LMM\Qc{p,s))m-p,s)) . (2.1) 



In the above expression Ei = Jmf + p^, M is the mock-meson mass, and the color 



wave function (subscript c denotes color) is written explicitly. Also, 0nLMi,(p) is a 
normalized momentum wave function, where n denotes all other quantum numbers of 
a state not connected to angular momentum (e.g., radial quantum number). The factor 



Y E appears due to our normalization convention for creation and annihilation 
operators |2^, {ha{k),h^^,{k')} = {2tt)^ ^ 6^(k — k')6aa', etc. The mock-meson states 
as given in (p.l|) are normalized to 2M. 



As already observed in [1^, the mock-meson approach suffers from a number of 
ambiguities, such as the choice for quark masses, or the definition of the mock-meson 
mass M. In the spirit of the method, the mock-meson mass should be defined as 



M = {Eg) + {Eq). However, as pointed out in |jT5[, the mock-meson mass has been 
introduced to give the correct relativistic normalization of the meson's wave function, 
and hence the use of the physical meson mass M instead of (Eg) + (Eq) may be more 
appropriate. We adopt the same approach, and write M = M. We also note that the 
heavier the mesons are, the less important it is how the mock-meson mass is defined, 
since the relativistic effects and binding energies become less significant. As far as quark 
masses are concerned, the self-consistency of the model requires the use of constituent 
quark masses. In our error estimates we have included variations of constituent light 
quark masses over a range of about 200 MeV, and also of heavy quark masses over 
a range of about 400 MeV, so that we believe that uncertainties introduced by a 
particular choice of quark masses are being properly taken into account. 

III. DECAY CONSTANTS 

Decay constants of heavy-light mesons are defined through matrix elements of vector 
and pseudo-vector currents between a meson state and the vacuum. Following 



standard definitions in the literature |Tl|], for pseudo-scalar (P), vector (V), scalar (S), 
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and pseudo-vector (A) mesons, we write 

(0|A^(0)|0-/2(A:)) = ^^^M^ (3.1) 

{0\V^mi-^,ie,k)) = ^_L_/^^^^Mv.,/,e'^ , (3.2) 

{0\V^{0)\Ot/,{k)) = ^^fsk^ , (3.3) 

(0|A'^(0)|l+2(e,/.)) = ^_L_/^^^^Ma^/,6^ , (3.4) 

{OmO)\lU^,k)) = j^jA,^.MA,^,e^ , (3.5) 

(0|^'^(0)|l3/2(6,fc)) = j7^Jv,;,My^,,e^ . (3.6) 



Note that in the heavy quark hmit ^"^^ ^1/2 states are degenerate (5- waves), and 
so are l^yg (-P-wave states). The spin 2 members of P-wave (13^2523/2) and 

D-wave (13^2,23^2) doublets do not couple leptonically due to conservation of angular 
momentum. 

In order to obtain expressions for decay constants in terms of integrals over 
momentum-space meson wave functions, we evaluate the matrix elements ( p.l| )-( pl6D 
in the meson rest frame using ( p.ll) . Of course, any choice of polarization for spin 1 
mesons should yield the same result. As mentioned earlier, current matrix elements 
between states defined in (|2.1| ) and the vacuum can be evaluated exactly with full Dirac 
spinors. Because of spherical symmetry, the momentum-space wave function can be 
written in the form 

(pnLMM = RnL{p)YLML{p) ■ (3.7) 

In the above Ylm^ are the usual spherical harmonics, and Rnbip) is the radial part of 
the wave function, where p denotes |p| henceforth. Using (|3.71 ), and keeping track of 
the relevant CG coefficients, we find that all heavy-light meson decay constants in the 
mock-meson approach can be written in the form 
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M 



(27r)3/2^ 



(3i 



where 



Fp{p) 



Fs{p) 



Fv,/M 



1 - 



1 + i 



mq + Eq){mQ + Eq) _ 

2 



p- 



1 1 



1 



("^Q + Fq) 
1 1 



2^2 1 



^no(p) , 

Rno{p) 
pRnlip) , 
pRnl{p) 



3 + 

2v^ 1 



pRnl{p) 
1 



P^Rn2{p) ■ 



Expressions (3J) and (|3.10|) were found in |12| and [17|, respectively. 

It is interesting to observe that in the hmit itlq oo (|3.8| )- (|3.14 ) become 

p^dp 



(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 



HL 



2^3 r— roo 

'Att 



(27r)3/2\ 



[mq + E g J 
2E, 



Ffip) 



(3.15) 



with 



F^'ip) 




= Rno{p) , 


FF'ip) 


= Fa,,Ap) 


= .pRnl{p) 
[mq + Eq) 




= , 






= . 





(3.16) 
(3.17) 
(3.18) 
(3.19) 



Equahty of fp and ^"^^ ^^^^ ^^^^ /s' ^^"^ /^i/2' ^^^^ vanishing of fA 



1/2 ' 



^13/2 



and /vg^j , are in the heavy quark hmit expected from the heavy quark symmetry. 
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IV. RELATIVISTIC QUARK MODEL 



In order to obtain numerical estimates for the decay constants of heavy-light mesons, 
we consider the simplest and widely used generalization of the non-relativistic Schro- 
dinger equation |^T],18-20| with Hamiltonian given by 



H 



+ p2 + + p2 + v{r) , (4.1) 



where for V{r) we take the QCD-motivated Coulomb-plus-linear potential 

4 (T 

V{r) = — — + br + c. (4.2) 
3 r 

For the sake of simplicity,^ we take as to be a fixed effective short range coupling 
constant. The effective string tension of the model can be determined from the re- 
quirement that the linear Regge structure of the model in the light-light limit agrees 



with the observed slope of the p trajectory p3|. Fixing m^ ^, other parameters can 
be chosen so that the model reproduces the observed spin-averaged spectrum of the 
known heavy-light states. One such set of parameters includes constituent quark masses 
mu,d = 0.300 GeV, m, = 0.483 GeV, = 1.671 GeV, and rub = 5.121 GeV, and also 
as = 0.498, b = 0.142 GeV^, and c = -0.350 GeV!^ As can be seen in Table ^ these 
parameters yield an excellent description of the observed spin-averaged heavy-light 
spectrum. 

We now turn to the discussion of pseudo-scalar and vector decay constants. Re- 
cently, Ref. [|l7i used (|4l|) with six different potentials, and with current quark masses 



^The running coupling constant was used in |11|. 

^This particular parameter set corresponds to the spin-averaged mass of the unknown Dq 



and D[ mesons (0^2 ^^"^ ^1/2 states) of about 2400 MeV. 
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from [El] , minimized the Hamiltonian with respect to the variational parameter (3 of a 



single harmonic oscillator (HO) wave function, 

and then used the wave function obtained in this way to get pseudo-scalar and vector 
decay constants from (|3.8|) , (p.9|) and (|3.10|) . However, a single harmonic oscillator 
(HO) basis state is not a suitable approximation for the meson wave function. Namely, 



lattice simulations pSf show that heavy-light wave functions fall exponentially with 
large r (~ e~^^'), and therefore HO wave functions (~ e"^^*"^/^) cannot be expected 
to reflect the correct dynamics of heavy-light mesons. If single basis states are used. 



a much better choice would be pseudo-Coulombic (PC) basis states p6| which fall 
exponentially with large r and appear to be in a good agreement with the lattice data, 
as can be seen in Figure |l]. 

Models such as the one we are using here are usually solved by diagonalizing the 
Hamiltonian matrix in a particular (truncated) basis, with basis states depending on 
some variational parameter As one increases the number of basis states, the de- 
pendence of eigenvalues and eigenfunctions on the variational parameter should vanish 
for the lowest states. In the case of QCD-motivated potentials the solutions obtained 
with the PC wave functions converge much more rapidly with an increase in the num- 
ber of basis states, than those obtained with the HO wave functions. We illustrate 
that in Figure 0, by plotting the dependence of energy of the lowest IS state on the 
variational parameter for = 1, 5 and 15 basis states, for both PC and HO wave 
functions. One can clearly see that the lowest IS* HO wave function is not a very good 
trial wave function in a variational calculation of ( |4. 1| ) (with QCD-motivated poten- 
tials). Furthermore, even if one believes that the = 1 HO result for a state energy 
is acceptable (it is roughly 50 MeV higher than the exact solution, as can be seen in 
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Figure H), that still does not justify the use of a single HO basis state as a meson wave 
function. This issue is clearly important in calculations where a correct description of 
meson dynamics is needed, such as calculations of meson decay constants. Results ob- 
tained by varying a single HO basis state are thus to be interpreted as non-relativistic 
estimates of some effective harmonic oscillator potential, and not as the results of a 
QCD-motivated relativistic quark model. 

One can now observe that if one uses enough basis states, the choice of basis wave 
functions should not matter, and pseudo-scalar and vector decay constants should be 
obtainable from the relativistic quark model considered here. The problem is, however. 



that the IS* wave function is divergent at the spatial origin p8|], i.e.. 

The singularity for r — > is related to the singularity of the short-range Coulomb 
potential. By increasing the number of (usually finite at r = 0) basis states, one 



is gradually beginning to see that singularity |2^. Furthermore, from ( [4. 41 ) one can 
see that the degree of divergence highly depends on the choice of a^. Because of 
that, one can expect that pseudo-scalar and vector decay constants cannot be reliably 
estimated within the model we are considering. In Figures ^ and |^ we demonstrate the 
dependence of the pseudo-scalar (D-meson) and vector (D*-meson) decay constants 
on the number of basis states (A^), for both PC and HO wave functions. As one can 
see, for small both fp and fvi/2 significantly increasing with an increase in 
A^. By including enough basis states, the dependence on A^ would eventually vanish.^ 



Because of the minus sign in ( |3.9D the results for fp are better behaved than those for fvi/2- 
For example, fp obtained with = 50 PC states are usually larger than those obtained with 
A^ = 25 by only a few MeV. On the other hand, the same increase in A^ in general leads to 
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However, as implied by ( |4.4| ), both fp and fv^^^ are quite sensitive to the particular 
choice of parameters of the model. In our calculations we have observed that results 
obtained with fixed can vary up to a few hundred MeV. Because of that, we were 
not able to obtain reliable estimates of fp and fvi/2 from the model considered in this 
paper .0 

One possible solution of the problem discussed above would be to replace the 1/r 
potential with the one-loop single gluon exchange potential, i.e., as — > (ys{r). The IS" 
solution of (^4.1|) in that case is still divergent, but the divergence is only logarithmic 
||28| . This should lead to much more stable results than the ones shown in Figures |^ 
and Ij. These results should also be much less dependent on the specific choice of the 
model parameters. In fact, such a calculation for fp (for D, Ds, B, and Bg mesons) 
was already performed by Capstick and Godfrey in [jl5| using the model of [|ll[]. The 
dependence of their results on the number of basis states was not shown, but the 



authors of |[T^ stated that they believed that the model overestimates pseudo-scalar 
decay constants (e.g., for D meson they found fp = 301 MeV with uncertainty of 
20%). Even though it is important to investigate what really happens with both fp 
and fvj^/2 s^ch a model, we shall not consider it in the present paper. 

We next discuss the heavy-light P and D-wave decay constants. While we were not 
able to obtain reliable results from ( [4.1|) and ( [4.2|) for the S-waves, the situation for P 
and D-waves is completely different. In Figures |, |, |^ and |] we show the dependence 
on the number of basis states (A^), for scalar (S), two pseudo- vector {A1/2 and A3/2), 



increase in fvi/2 by several hundred MeV. 

^From Figures 1^ and 1^ it should be clear that in the model considered here the ratio fp/ fvi 
also cannot be determined with reasonable errors. 
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and vector (V3/2) decay constants, respectively. All the results shown are for the D** 
mesons. As one can see in those figures, in general only a few basis states are needed 
for results to become independent of A^, even though the derivatives of the actual IP 
and ID wave functions are singular at spatial origin I^Oj-Q Furthermore, as A^ increases 
the HO results approach the PC results (always from below) which shows that the 
difference between the two basis sets is slowly vanishing. However, even with 15 basis 
states (when the state energy obtained from the model is essentially equal for both 
PC and HO wave functions), we can still see the difference for (Figure |]) and for 
(Figure 0). This reflects the difference in the wave functions obtained from the 
two basis sets. The reason why both PC and HO basis states yield almost the same 
results for fs (Figure |]), even though 0^^2 state is also a P-wave, is the minus sign in 
( p.ll| ). Of course, because of the much more rapid convergence, the PC results are to 
be preferred over the HO results. 

Our calculations of P and D-wave decay constants showed that their dependence 
on the particular choice of the model parameters is signiflcantly smaller than the cor- 
responding dependence of fs and /v^/j- present the results for D** , D**, B**, 
and B** mesons in Tables Q, [T], |rV|, and |V], respect ively.[| To obt ain these results 
the effective string tension b of the model was determined from the observed slope of 



"By fixing all input parameters, the sensitivity of the decay constants on the number of 
basis states was investigated. To achieve an accuracy of 0.1 MeV for fs and /vg^j ^ fittle 
as 10 PC basis states usually were needed, while to achieve the same accuracy for fAi/2 ^-i^d 
requires in general about 50 to 75 PC basis states. 

'^All results given in Tables |l| through were obtained with 25 PC basis states, which was 
more than enough for the accuracy of less than 1 MeV in all cases considered. 
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the p trajectory. For a fixed 171^4 other parameters were obtained from the spectrum 
of the known heavy-light states. Experimental meson masses were used in ( p.8[ ) only 
when their quantum numbers were unambiguously determined. Else, we used model 
predictions for the appropriate spin-averaged masses, which are also shown in Tables 
through |V[ 

In order to estimate uncertainties introduced by a particular choice of the con- 
stituent mass of u and d quarks, we have varied 171^4 in the range from 150 MeV to 
350 MeV . For a given 171^4-, by adjusting c we have also varied constituent heavy quark 
masses in the range of about 400 MeV (e.g., tjIc was varied in the range from about 
1.3 GeV to about 1.7 GeV). We emphasize that a good description of the spin-averaged 
heavy-light meson spectrum was always maintained. 

Results for the decay constants obtained in this way depend on the assumption for 
the unknown spin-averaged mass of and D[ mesons (0^^2 ^"^^ I1/2 states). To take 
into account ambiguities introduced in our results in that way, we have repeated all 
calculations for this unknown mass in the range from 2200 MeV to 2450 MeV . Errors 
quoted in Tables |I| through |V] reflect the uncertainty due to the unknown P-wave mass, 
as well as the uncertainties related to the choice of constituent quark masses discussed 
above. 

As one can see from those tables, in spite of the fact that our calculations are 
performed for a broad range of constituent quark masses, and also for a wide range of 
the unknown P-wave mass, as long as a good description of the observed heavy-light 
meson spectrum is maintained, the P and D-wave heavy-light decay constants are all 
predicted rather precisely. It is also interesting to observe that the decay constants of 
strange 0j*y2) I1/25 Ig^g states are slightly smaller than those of the corresponding 
non-strange states. The main reason is (besides the 1/\/M dependence of ( |3.8| )) the 
light quark dependence of ( |3.11D , (|3.12| ) and (|3.14| ). On the other hand, ( p. 131) does 
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not depend on the light quark mass, so that ^{mq + Eq)/Eq factor in ( |3.8D plays a 
much more significant role, and as a result /A3/2 ^^e strange states are larger than 
the ones for non-strange states. Also note that fs for and Bsq are larger than those 
of the corresponding Dq and D^o states, while it is the other way around in the case of 
"^^^ reason for this are the minus and plus signs in and (|3.12| ), respectively. 

Finally, the fact that I'^j^ and l^j^ B** states have decay constants smaller than those of 
the corresponding D** states, can be easily explained with the l/(mQ+i?g) dependence 
of ( CT) and ( CT )- 

V. CONCLUSION 

In this paper we have examined decay constants of heavy-light mesons within the 
mock-meson approach ||7|-pil[. We obtained all the relevant expressions in the j-j 



coupling scheme. For numerical estimates we employed a simple and widely used 
relativistic quark model It is based on a spinless Salpeter equation with 



QCD-motivated Coulomb-plus-linear potential. The effective string tension is chosen 
so that the Regge structure of the model in the light-light limit is consistent with 
experiment, and other parameters are based on the good description of the known 
spin-averaged heavy-light meson masses. 



Due to the singular nature of the L = wave functions at spatial origin [2S], we were 
not able to obtain reliable estimates of pseudo-scalar and vector decay constants. On 
the other hand, even though we have allowed for large variations of input parameters, 
our results show that the model predicts a rather narrow range for all lowest P and 
D-wave heavy-light decay constants. 

Such precisely predicted decay constants allow us to estimate the D**^ production 
fractions in b decays governed by the b ccs' transitions under the factorization 



14 



assumption. Quantitative predictions regarding the interference of color-allowed and 
color-suppressed amplitudes in —>■ D**^ {n^ , , ai , . . .} modes can now be formu- 
lated. These and some other consequences of our findings will be discussed elsewhere 
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TABLES 

TABLE I. Relativistic quark model predictions compared to experimental spin-averaged 

heavy-light meson masses. Parameters of the model are m„ = 0.300 GeV, nis = 0.483 GeV, 
rUc = 1.671 GeV, nn, = 5.121 GeV, as = 0.498, b = 0.142 GeV"^, and c = -0.350 GeV. The 
unknown Dq and D[ mesons {0^2 ^^'^ ^1/2 states) were assumed to have a spin-averaged 
mass of 2400 MeV. Heavy quark symmetry arguments then lead to the spin- averaged mass 
of 2502 MeV for the corresponding Dgo and -D^^ mesons. 



Meson 


State 




Experiment 

[MeV] 


Theory 

[MeV] 


Error 

[MeV] 


£>(1867) 


0^/2 


> 


1S'(1974) 


1971 


-3 


D*(200^) 


^1/2 










Dn(r^ 2400) 


Own 
"1/2 










£>i(~ 2400) 


1 + 

1/2 


> 


1P(2431) 


2434 


+3 


D^ f2425'l 


1+ 










L»|(2459) 


2+ 

3/2 










£>s(1969) 


1/2 


> 


15(2076) 


2079 


+3 


r>*(2112) 


^1/2 










D,o(- 2502) 


0^-/2 










i:>^i(~ 2502) 


1 + 
^1/2 


> 


1P(2540) 


2537 


-3 


£)si(2535) 


1 + 
^3/2 










Z)*2(2573) 


2+ 

^3/2 ^ 










5(5279) 


0V2 


> 


15(5314) 


5314 


+0 


5* (5325) 


^1/2 










5s (5374) 


Ov2 


> 


15(5409) 


5409 


-0 


5* (5421) 


11/2 ^ 
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TABLE II. Decay constants of heavy-light D** states, as obtained from the relativistic 
quark model. Whenever possible we used experimental meson masses. If these were unknown, 
we used model predictions for the spin-averaged masses. 



Meson 


State 


fi 

[MeV] 


Dfl867^ 

i:)*(2009) 

Do (2410 ± 40) 

Di(2410±40) 

Di(2425) 

Di'(2700 ± 55) 


IS 07/0 

' 1/2 
IS, 1^/2 
IP, 0+2 

IP 1 + 

^3/2 

ID, I3/2 


not reliable 
not reliable 
139 ± 30 
251 ± 37 
77± 18 
48 ±7 


TABLE III. Decay constants of heavy-light D** states, as obtained from the relativistic 
quark model. Whenever possible we used experimental meson masses. If these were unknown, 
we used model predictions for the spin-averaged masses. 


Meson 


State 


fi 
[MeV] 


D5(1969) 

D*(2112) 

D,o(2510 ± 45) 

D^i(2510±45) 

Dsi(2535) 

D^'i(2795±55) 


15, 0-/2 

IS, 1^/2 

IP, 0+2 
1-^' 13/2 

ID, 13/2 


not reliable 
not reliable 
110 ± 18 
233 ± 31 
87± 19 
45 ±6 
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TABLE IV. Decay constants of heavy-light B** states, as obtained from the relativistic 
quark model. Whenever possible we used experimental meson masses. If these were unknown, 
we used model predictions for the spin-averaged masses. 



Meson 


State 


[MeF] 


^('5279') 

S*(5325) 

So (5765 ± 60) 

5i(5765±60) 

Si(5765±60) 

Si'(6040±70) 


IS or,o 

' 1/2 
IS, 1^/2 
IP, 0+2 

lV2 
^3/2 

ID, 1-/2 


not reliable 
not reliable 
162 ± 24 
206 ± 29 
32± 10 
18 ±3 


TABLE V. Decay constants of heavy-light B** states, as obtained from the relativistic 
quark model. Whenever possible we used experimental meson masses. If these were unknown, 
we used model predictions for the spin-averaged masses. 


Meson 


State 


fi 
[MeV] 


5,(5374) 

S*(5421) 

5,0 (5860 ± 65) 

5^i(5860±65) 

5,i(5860±65) 

5^\(6130±75) 


IS, 0-/2 

IS, 1^/2 
IP, 0+2 

1^' lV2 

IP 1 + 

^3/2 

ID, 1-/2 


not reliable 
not reliable 
146 ± 19 
196 ± 26 
36± 10 
17 ±3 
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FIGURES 

FIG. 1. Comparison of the pseudo-Coulombic (PC, Ris{r) = 2/?^/^e~^'"), and the har- 
monic oscniator (HO, Ris{r) = 2/33/2/vri/4e-/32rV2)^ ig configuration space wave functions 
with the lattice data [^5|. For both PC and HO wave functions we used /? = 0.40 GeV. 

FIG. 2. Convergence of the 15 state mass of ( p?TD and (^^21), with mi = mu,d = 0.300 GeV, 
m2 = rric = 1.671 GeV, b = 0.142 GeV^, c = -0.350 GeV, and a, = 0.498. 
Pseudo-Coulombic (PC, full lines) and harmonic oscillator (HO, dashed lines) wave func- 
tions with = 1, 5, and 15 basis states. 

FIG. 3. Dependence of the pseudo-scalar {D meson, 0~^^ state) decay constant fp on the 
number (A^) of pseudo-Coulombic (PC), and harmonic oscillator (HO) basis states. We have 
used parameters given in the text, i.e., mi = ruu^d = 0.300 GeV, in2 = rric = 1.671 GeV, 
b = 0.142 GeV'^, c = -0.350 GeV, and a, = 0.498. 

FIG. 4. Dependence of the vector {D* meson, l~y2 state) decay constant fvi/2 
number (N) of pseudo-Coulombic (PC), and harmonic oscillator (HO) basis states. We have 
used mi = m^^a = 0.300 GeV, m-i = m^ = 1.671 GeV, b = 0.142 GeV^, c = -0.350 GeV, 
and as = 0.498. 

FIG. 5. Dependence of the scalar {Dq meson, Oy^ state) decay constant fs on the number 
(A^) of pseudo-Coulombic (PC), and harmonic oscillator (HO) basis states. We have used 
mi = mu4 = 0.300 GeV, m2 = m^ = 1.671 GeV, b = 0.142 GeV"^ , c = -0.350 GeV, and 
as = 0.498. 
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FIG. 6. Dependence of the pseudo-vector {D[ meson, ly^ state) decay constant /ai/2 
on the number (N) of pseudo-Coulombic (PC), and harmonic oscillator (HO) wave func- 
tions. We have used mi = mu,d = 0.300 GeV, m-i = = 1.671 GeV , b = 0.142 GeV'^, 
c = -0.350 GeV, and = 0.498. 

FIG. 7. Dependence of the pseudo-vector {Di meson, 1'^^^ state) decay constant /A3/2 
on the number (N) of pseudo-Coulombic (PC), and harmonic oscillator (HO) wave func- 
tions. We have used mi = mu,d = 0.300 GeV, m2 = rric = 1.671 GeV, b = 0.142 GeV'^, 
c = -0.350 GeV, and = 0.498. 

FIG. 8. Dependence of the vector {D" meson, 13^2 state) decay constant fv^/^ '^^ the num- 
ber (N) of pseudo-Coulombic (PC), and harmonic oscillator (HO) wave functions. We have 
used mi = m„,d = 0.300 GeV, m2 = rric = 1.671 GeV, b = 0.142 GeV"^, c = -0.350 GeV, 
and as = 0.498. 
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